Many-body treatment of quantum transport through single molecules 
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We develop a theoretical approach to quantum transport through single conjugated organic 
qq ■ molecules that accurately accounts for all molecular excited and ground states via a semi-empirical 

f*^ many-body molecular Hamiltonian. We then calculate the linear and non-linear transport proper- 

ties of a benzenedithiol molecule covalently bonded to two gold electrodes at room temperature, 
£Sj ■ eliminating all free parameters by comparing the calculated thermopower and linear-conductance 

' with experiment. The resulting 'molecular diamond' structure exhibits experimentally observed 

' many-body effects inaccessible to mean-field approaches. 

Q: 

' Electron transport through single molecules is of interest both fundamentally and as a means to probe numerous 
y—{ \ important chemical 1 and biological processes^, with myriad potential device applications 3 . In spite of the growing body 
of experimental evidence indicating that intra-molecular interactions play an important role in transport problems 
1 ^ '' involving single-molecule junctions^£&£, theoretical studies of these junctions typically focus on either many-body 
O . phenomena, such as the Kondo effect^, while neglecting the internal structure of the molecule, or strive for a realistic 
description of the molecule's electronic structure via ab initio methods 3 -^, which treat interactions at the mean-field 
, level. 

In this article, we develop a many-body non-equilibrium Green's function (NEGF) approach to the problem of quan- 
tum transport in multi-terminal single-molecule devices. The electronic degrees of freedom most relevant for trans- 
■ port, i.e., the 7r-electrons in conjugated organic molecules, are modeled via a semi-empirical Pariser-Parr-Pople (PPP) 
Hamiltonia n 10 ' 11 ! 12 . The molecular Hamiltonian is solved via exact diagonalization, thereby including intra-molecular 
correlations non-perturbatively. The coupling of the molecule to the macroscopic leads is treated perturbatively within 
the NEGF formalismi 3 .^^. 

Electronic transport is a property of a junction: molecule plus leads. The Hamiltonian of a junction may be written 
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as H — H mn \ + -ffip ar i s + Htnnt where the first term is the 7r-electron molecular Hamiltonian^ 
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, where S na creates an electron of spin a in the 7r-orbital of the nth carbon atom, e„ is the on-site orbital energy, t nm are 
the tight-binding matrix elements between orbitals n and to, U nm are the intra- and intersite Coulomb interactions, 
, and ( ) indicates a sum over nearest- neighbors. 

The effective charge operator 1 - for orbital n is given by 
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. where C' na is the capacitive coupling between orbital n and lead a, e is the electron charge, and V a is the voltage 
on lead a. The lead-molecule capacitance values are elements of a full capacitance matrix^ C , which includes 
both the intra-molecular and inter-molecular capacitances. This full capacitance matrix obeys the usual zero-sum 
rules required for gauge invariance and is constructed in accordance with the intra-molecular interaction energies (i.e. 
C™ = e 2 U~m) and the geometry of the junction, where we make the ansatz that C na is inversely proportional to the 
lead-site distance. 

The two remaining terms in the junction's Hamiltonian, -ffi ea ds and H tU n, account for the isolated leads and lead- 
molecule coupling, respectively. Due to their large density of states, the macroscopic metallic leads (labeled a 6 
[1, . . . , M]) may be modeled as non-interacting Fermi gases, with tunneling matrix elements V n k from a level fc of 
energy within lead a to a nearby 7r-orbital n of the moleculeJ^ These terms are given explicitly by 
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where ct is a lead-electron creation operator. 

The effective interaction energies for a 7r-conjugated polymer are given b y 10 ' 11 ' 12 
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where U is the on-site Coulomb energy, e is the dielectric constant, a = (t//14.397eV) , and R nm is the distance 
between orbitals n and m. The dielectric constant e accounts for screening effects due to both the c-electrons and 
any environmental considerations, such as any non-evaporated solvent, 12 The solutions of the PPP Hamiltonian, used 
in conjunction with these experimentally derived interaction energies, accurately account for all ground and excited 
state energies and wavefunctions of the free molecule^ 2 -. 

When the isolated molecule, with its discrete molecular energy spectrum, is coupled to macroscopic metallic leads 
a junction is formed and the bare molecular energy levels become dressed by interactions with the metal's electronic 
state continuum. These interactions give rise to width broadened and shifted molecular energy levels which are 
imbued, correspondingly, with finite lifetimes. We are interested in the non-equilibrium steady state behavior of these 
junctions both as a function of bias voltage and of molecular gating voltage. 

A considerable simplification to the problem of quantum transport in this fundamentally open system is obtained 
through the use of the Green's functions G(E) and G < (E), which are Fourier transforms of the retarded and Keldysh 
Green's functions, G W(W (t) = -*0(t)({d»«r(t)>4»o*(°)}> and G< lrTmrT ,{t) = i{d\ nrT , (0) d na {t)), respectively! 3 .^. 
The intuitive interpretation of the retarded Green's function as the total probability amplitude for a system with a 
particle (electron or hole) at site n to be found a time t later in a state with an added particle at site m, is its principle 
advantage in solving quantum transport problems. For an isolated molecule, we derive the (retarded) Green's function 
from the molecular Hamiltonian and find 
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where \v) and \v') are many-body eigenstates satisfying H nlo \\v) — E v \v), and V(y) is the probability that the isolated 
molecule is in state \v). In linear response, the probabilities V(v) in Eq. ^ are determined by equilibrium thermo- 
dynamics, Viv) = e~P( El '~ flN " S> /Z, where Z is the grand partition function of the molecule at inverse temperature /3 
and chemical potential /i. For nonlinear transport, V(y) is given by the appropriate limit as V n k — > 0. This is the 
limit of sequential tunneling, where V(y) can be determined from detailed balance in a rate-equation approach 19 . 
The Green's functions of the full system are 
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where £ and S < are the self-energies describing the dressing of the bare molecular states by the lead-molecule coupling. 
Given the Green's functions, the current flowing into lead a is then 13 



J dE T a (E) lm[G<^ aa (E) + 2f a (E)G aa , aa (E)}, 



(8) 



where T a (E) = 2nJ2kea \Vnk\ 2 (>(E - efe) is the tunneling width for lead a, f a {E) = {1 + cxp[(£' - ^ Q )/fc B T]} _1 , and 
a is the 7r-orbital coupled to lead a. 

It is nontrivial to derive exact expressions for the self-energies describing lead-molecule coupling. Wick's theorem 
cannot be used in the usual way to derive perturbative corrections because -ff mo i is not quadratic in fermion operators. 
Instead, diagrammatic expansions must be done in the independent-particle basis, where the rightmost term in Eq. 
(fTj) is also treated as a perturbation. To lowest (zeroth) order in U nm , the self-energies are simply 
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Even in this lowest-order expression for the self-energy, the lead-molecule coupling is still treated at the same level 
of approximation as in mean-field ah initio models^, where lead-molecule correlations are neglected. Eqs. ([ MTU]) 
represent a conserving approximation; current is conserved and the spectral function obeys the usual sum-rule. 
Transport is treated at the level of elastic cotunneling theory- 2 ^, while sequential tunneling results^ are recovered 
in the limit T a /T — > 0. This approximation is valid when there are no unpaired electrons in the molecule or at 
temperatures large compared with the Kondo temperature.— 

Within this approximation, we are left with a key result: G mo i alone encapsulates the intra-molecular correlations. 
The full NEGF current expression (jSj) then reduces to the well-known multi-terminal Buttiker formula 21 
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is the transmission probability from lead a to lead (3 and a (b) is the orbital coupled to lead a (/?) . 

Given the large uncertainties concerning the lead-molecule coupling 22 we shall take the so-called broad-band limit 14 
in which the T a are taken as constants characterizing the lead-molecule coupling. Typical estimate s 23 ' 24 ' 25 using the 
method of Ref. 26 yield r Q < leV for organic molecules coupled to gold contacts via thiol groups. In the broad-band 
limit, we find that the non-equilibrium weighting function becomes V [y) = J2a=i ^aPa i. v ) I T Q , where V a [y) 

is the equilibrium probability for a state v in the grand-canonical ensemble at the temperature and chemical potential 
of lead a. 

As a first application of our approach, we consider the well studied two-terminal 1,4-benzenedithiol (BDT) molec- 
ular junction with gold leads i 27 ' 28 ' 29 ' 30 ' 31 Systematic molecular energy-level studies of experimental condensed-phase 
benzene spectra have found the parameterization t nm — 2.64eV, U = 8.9eV, and e = 1.28 of Eqs. (P) and ([5]) to 
give the best quantitative agreement for all ground and excited states^ 2 Simulations discussed in this paper were all 
performed with symmetric capacitive and lead- molecule couplings for a junction at room temperature (T=273K). 

The transmission probability given in Eq. (|12p is also the linear-response conductance in units of the conductance 
quantum Go = 2e 2 /h. This conductance is a function of both the lead-molecule tunneling width, T, and the difference 
between the lead's Fermi level and the isolated molecule's HOMO-LUMO gap center, [iq. Our model therefore has 
two phenomenological parameters which we shall determine by fitting to experiment. 

Thermopower (Seebeck coefficient) measurements 32 offer an experimentally accessible means to determine 3 - 3 - the 
lead-molecule chemical potential mismatch, /ia u — Mo, where /iq = (chomo + ^lumo) /2 and chomo {^lumo) rep- 
resents the HOMO (LUMO) energy level. Although usually associated with bulk materials, the Seebeck coefficient, 
S, can also be obtained for a molecular junction by measuring the voltage difference created across the junction in 
response to a temperature differential. Beginning from Eq. ((11)) and the definition of thermopower as AV = — SAT, 
where V is the voltage and T is the temperature, one finds 34 

where fj, = /j, a = fip is the chemical potential of both leads. An important feature of the Seebeck coefficient is that 
asymmetries in conductance with respect to /ja u — fJ-o manifest themselves in the sign and magnitude of the above 
integral, thereby lending insight into the nature of the charge transport in an energy barrier dependent, but molecule 
number independent manner. 32 ' 33 Reddy et al^ have measured the Seebeck coefficient of a BDT-junction and found 
it to be (8.7 ± 2.1) (xV/K. By equating this experimental value with that calculated using Eq. (fl3|) we find, as shown 
in Figure [TJ that ^Au — Mo = — (3.5±0.25)eV. This value indicates that the Fermi level of the gold lies about 1.5eV 
above the HOMO resonance, which validates the notion that these devices are dominated by p-type (hole) transport. 

The only other free parameter is the lead-molecule coupling, T, which can be found by matching the linear-response 
conductance to experiment. Although there is a large range of experimental values^ the most reproducible and 
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FIG. 1: The Seebeck coefficient as a function of the lead-molecule chemical potential mismatch for BDT-junction. Inset: A 
closer look at the curve showing, away from a resonance, that the Seebeck coefficient is only weakly dependent on the lead- 
molecule coupling F. The horizontal line represents the mean experimental value of .— We find that fj.Au is 3.5eV 
below the center of the LUMO-HOMO gap (/x ) of BDT 

lowest resistance contacts were obtained by Xiao, Xu and Tao2^, who report a single-molecule conductance value of 
O.OllGo- This fixes r=(0.53±0.11)eV, within the range predicted by other groups for similar molecules. With the 
final parameter fixed, we can now use our many-body approach to predict the linear and nonlinear response of this 
molecular junction. 

Initially, a calculation of the transmission probability T{E)\e =1 i as a function of /x— /itAu was performed for the BDT- 
junction by exact diagonalization of Eq. JT]) and the construction of the requisite Green's functions. The chemical 
potential /x is related to the junction's gate voltage V g via the capacitance C g . The resulting electron addition 
spectrum, shown in the top part of Figure [21 exhibits several striking features arising from our many-body treatment: 
an irregular Coulomb blockade-like peak spacing, indicative of charge quantization; an increased LUMO-HOMO gap, 
reflecting incipient Mott-Hubbard correlations; and non-symmetric, Fano-like resonance lineshape arising from the 
non-zero thermodynamic weighting for multiple states near a resonance. Further, since the thermal energy fc^T is far 
less than any molecular electronic resonance or T, the zero bias addition spectrum is also a map of the ground state 
electronic structure of BDT. Our approach therefore reproduces the key results of both the Coulomb blockade and 
resonant tunneling regimes simultaneously, including intramolecular many-body correlations exactly. 

We next calculate the differential conductance, dl/dVbiaa, as a function of f\i and Vbi as - The results of this calculation, 
shown in the lower part of Figure [U exhibit a 'molecular diamond' structure which is a molecular analog of the 
Coulomb-diamond measured in many quantum dot transport experiments. 3 — Apart from the central LUMO-HOMO 
gap, the resonance spacing in the top half of Figure [5] can be roughly explained via a capacitive model in which the 
molecule is characterized by a single capacitance C = e 2 / (U nm ), where () indicates an average over site number. With 
the aforementioned parameterization 12 - (U nm ) ~5.11eV whereas the LUMO-HOMO gap ~10eV, indicating significant 
deviations from a simple capacitive charging model due to molecular structure (note also the enlarged diamonds with 
charge ±4). While charge quantization has not yet been observed in BDT junctions due to the difficulty of gating such 
small molecules 3 -^, the unambiguous observation of Coulomb blockade in junctions involving larger molecules^&i' 3 -^, 
with smaller charging energies, indicates that such interaction effects, lying outside the scope of mean-field approaches, 
are undoubtedly even more important in small molecules, like benzene. 

A cross section at zero bias reproduces the addition spectrum shown in the top part of the same figure. Increasing the 
bias voltage, we find that the resonances split into negatively sloped particle- like (\e\dV/dfi = — {C\ + C2 + C g ) jC\) 
lines and positively sloped hole-like (\e\dV/dfi = + (Ci + C2 + C g ) /C2) lines, where Ci and C2 are the mean lead- 
molecule capacitances. In the symmetric coupling case (Ci = C2 with C g jC\ <C 1) the lines therefore have slopes of 
-2 and +2 for particle-like and hole-like lines, respectively. 

Within the V-shaped outline traced by the particle- like and hole- like lines we find evidence for excited state transport 
in the many narrow, nearly parallel resonance lines. The structure of these resonances is highly dependent on particle 
number, as one would expect from the underlying many-body symmetries. Excited state transport has been observed in 
quantum dot 3 -^ and nanotube experiments 3 -^ but hasn't yet been unambiguously identified in single-molecule junctions. 

In conclusion, our theoretical approach combining a semi-empirical PPP Hamiltonian with NEGF methods provides 
a computationally tractable means to study intra-molccular many-body correlations and excited state transport in a 
quantitatively accurate way. Applying our theoretical approach to a BDT-junction with two gold leads, we solved for 
the only two free parameters in our model, finding /ia u — Mo = —(3.5 ± 0.25)eV and T = (0.53 ± 0.11)eV, indicating 
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FIG. 2: The linear and non-linear response of a BDT-junction with two gold leads (^Au — fJ-o = — 3.5eV). Top: The electronic 
addition spectrum as a function of gating potential at zero bias exhibits non-symmetric, Fano-like resonance lineshapes and a 
widened LUMO-HOMO gap due to Mott-Hubbard correlations. This is also the linear-response conductance in units of the 
conductance quantum Go = 2e 2 jh. Bottom: Each column is the logarithm of the differential conductance, log 10 [(dl/dV) /Go], 
calculated for a given gating potential fi — /xau- The number placed inside each constant-particle- number diamond represents 
the molecule's charge state, with respect to neutral benzene. The conductance resonance lines with positive slope indicate 
hole-like excitations while the negative slope indicate particle-like transport, which at zero bias are degenerate. The finer, 
nearly parallel lines inside the outline of the particle- like and hole-like resonances arise from excited state transport. 

that the junction's charge transport is dominated by holes (p-type transport) and that the lead-molecule coupling is 
within the expected range. At non-zero bias, the differential conductance of a BDT junction is predicted to exhibit a 
'molecular diamond' structure analogous to the 'Coulomb diamonds' observed in quantum dots^, but with significant 
modifications due to molecular structure. 

The authors acknowledge useful discussions with S. Mazumdar and D. M. Cardamone. 
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